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Abstract 

We investigate the perturbative (in gsNos) backreaction of localized D8 
branes in D4-D8 systems including in particular the Sakai Sugimoto model. 
We write down the explicit expressions of the backreacted metric, dilaton and 
RR form. We find that the backreaction remains small up to a radial value of 
u <S 4/(3s^D8)) and that the background functions are smooth except at the 
D8 sources. In this perturbative window, the original embedding remains a 
solution to the equations of motion. Furthermore, the fluctuations around the 
original embedding, describing scalar mesons, do not become tachyonic due 
to the backreaction in the perturbative regime. This is is due to a cancela- 
tion between the DBI and CS parts of the D8 brane action in the perturbed 
background. 
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1 Introduction 

By now, AdS/CFT has become a standard tool in theoretical physics for the study 
of gauge theories at strong coupling. In many "stringy" models of gauge dynamics 
fundamental matter is included by embedding a set of "flavor branes" in addition to 
the "glue/color branes." In such a setup, the strings connecting only to the "glue 
branes" are in the adjoint of the U(N C ) group, giving gauge particles (multiplets) , 
and those connecting only to the "flavor branes" are in the adjoint of U(Nf), giving 
the mesons (meson multiplets), and those connected to both the "flavor" and "glue" 
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branes are in the fundamental representation of both the U(N C ) and U(Nf) groups, 
giving the quarks (matter /quark multiplets). Anti-quarks are obviously depicted by 
similar strings with the opposite orientation. 

In principle for large N c and large Nf one could go to a combined near horizon 
limit, translate the branes into fluxes and derive the gravity background that is a 
holographic dual of a gauge system with gluons and quarks ( and, if the model is 
supersymmetric, their supersymmetric partners). However, in practice such models 
are being constructed using the probe approximation. In this approximation one 
uses a gravity background built from the near horizon limit of a large N c glue branes 
and adds to it a set of Nf flavor probe branes. The basic assumption of the probe 
approximation is that for the case of Nf <C N c the backreaction of the probes on the 
background can be neglected. The flavor physics is then extracted by analyzing the 
effective action that describes the flavor branes in the glue background, namely the 
DBI action plus the CS action. This practice was introduced in [1] in the context of 
the AdS^ x S 5 model, for a confining background in [2] and subsequently in a large 
variety of other models [3], [4], [5]. Exceptions to this are certain fully backreacted 
non critical models like [6], [7], [8], and [9]. 

A landmark holographic model of chiral symmetry and chiral symmetry breaking 
is the model of Sakai and Sugimoto [10]. This model is based on the incorporation of 
a stack of Nf D8 and Nf anti- D8 flavor branes into the background of near extremal 
DA branes [11]. In the latter background one compactifies one of the world volume 
coordinates of the D4 branes on a circle of radius R. For energies E -C 1/R the 
background describes a four dimensional system with gluon degrees of freedom plus 
contaminating Kaluza Klein modes. The profile of the flavor branes determined by 
the DBI action is that of a U shape. This provides a simple geometrical picture 
of chiral symmetry breaking, namely, for large radial direction (see figure [2]), which 
corresponds to the UV limit of the gauge theory, the stack of the D8 branes and of 
the anti-D8 branes are separated and hence there is a U^Nf) x Un(Nf) symmetry, 
and in the IR limit the two stacks merge one into the other and thus only the 
diagonal U(Nf) survives as a symmetry. A variety of physical properties of meson 
and baryon physics has been extracted from the model. These include the massive 
meson spectrum, the massless Goldstone pions [10], certain decay rates [12] as well 
as the thermal behavior of hadrons [13], [14]. 

The validity of the model [10], is the same as all other probe models: Nf <C N c . 

bearing in mind that one must keep the open string spectrum associated with the flavor branes, 
even though they are producing macroscopic flux. Such a case is a full backreaction, but not a 
decoupling limit. Recall that global symmetries of the field theory translate into gauge symmetries 
in the gravity: gauging a large Nf flavor group using supergravity alone is unfeasible. 
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To contact to real hadron physics, one obviously is interested in the case where the 
number of flavors is similar to that of the colors and both are not large. To get down 
to small N c one will have to invoke a full string theory rather than a effective gravity 
model. However, increasing the ratio of Nf/N c can still be done in the context 
of an effective field theory, provided we go beyond the probe approximation and 
incorporate the backreaction of the flavor branes on the gravity background. This 
may enable us to determine the flavor dependence of certain physical properties of 
the gauge theory which we expect to be Nf dependent, for example the beta function, 
or the ratio of viscosity to entropy density of the quark-gluon fluid [15] o 

Similar studies for localized backreactions in D3-D7 systems include [17-21]. In 
[17], very general framework for studying type IIB supergravity with metric/five form 
and holomorphic dilaton/axion. The work of [18-20] studied when the D7 branes 
are located at singular points in manifolds, and [21] studied the form of the solution 
to the equations following from the D3-D7 system and effects on other probe branes 
in such backgrounds. All these studies worked using the supergravity alone, while 
here we will derive delta function source terms from an action of the form 

SBuik + K 8 (Sdbi + Scs)- (1-1) 

We will use this action to determine how to source the bulk fields. Although we 
obtain the full equations of motion from this, we will study their solutions in a 
perturbative limit. Therefore, while we are going beyond the probe approximation, 
we are still confined to the regime Nf <C N c for the simple reason that we want the 
series in powers of Nf/N c to converge quickly. 

In fact there is an even more important motivation to explore the model of [10] 
beyond the probe approximation, and that is the issue of the stability of the model. 
One may wonder whether the model is stable only in the probe approximation and 
that the backreaction of the probe branes on the background does not destabilize 
the setup. A simplified picture of the model is that of the circle of the compactified 
direction with the two endpoints of the stacks of probe branes and anti branes which 
can be represented as a +Nf charge located at one point on the circle and — Nf 
charge located at the antipodal point. In this simplified "electrostatic" setup if one 
of the charges gets a slight perturbation in one direction it will be attracted to the 
opposite charge and will not be driven back to its original location. Moreover, the 
antipodal setup described in [10] has been generalized to a family of setups where 

2 In [16], the leading order correction in Nf/N c of this property was determined in the context 
of a model with D7 branes in near extremal Ads 5 x S 5 background. However, in this model it was 
shown that the zero mode, which is equivalent to smearing, is all that is necessary to this level in 
N f /N c . 
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the separation distance between the brane and anti-brane is taken to be L < irR. 
For these cases the "electrostatic instability" is even more severe. The question is 
therefore whether this naive intuition is justified and the backreaction of the probe 
branes indeed destabilizes the model. On the other hand there is a naive argument 
why the perturbative backreacted system should be stable and non tachyonic. Since 
the gauge holographic dual of the model before purturbing it has a spectrum with 
a mass gap ( apart from the pions), a small a small perturbation cannot bridge the 
gap and produce tachyonic modes [22]. 

Further, one may wonder what happens to the the dilaton tadpole condition, 
given that this is a D8 D8 on a circle, and both branes and anti branes source 
the dilaton in the same way. Hence, for these codimension one flavor branes one 
anticipates that the dilaton will have a cusp behavior at the location of the probe 
branes as well as a cusp (and not an anti-cusp) at the anti-brane. It seems naively 
that there is no way to sew together these two cusps. 

Thus the goal of this paper is to compute the leading order backreacted back- 
ground, and address the stability, and the dilaton tadpole. To do so, we write down 
the full action of the system (of the form in fll.lR ) which is composed of the action 
of massive type 11a supergravity and the nine dimensional DBI +CS actions asso- 
ciated with the the D8 flavor branes. At this point one often invokes a smearing 
of the flavor branes along their transverse direction [23-25] which renders the com- 
bined action into a ten dimensional one. This approach simplifies the analysis by 
turning the equations of motion (EOMs) into ordinary differential equations (ODEs) 
of some radial variable. However, we expect that certain physical questions may not 
be answered using this procedure, for example the stability discussed above. Thus 
we avoid using the smearing technique and we keep the flavor branes as localized 
objects. This yields delta function source terms for the equations of motion of the 
graviton, dilaton and the -F^io) RR field strength form associated with the D8 branes. 
The complexity of these equations is increased, relative to the smearing technique, 
because the EOMs are now partial differential equations (PDEs); the relevant func- 
tions must depend on the coordinate(s) transverse to the flavor brane. |f| 

We solve these equations perturbatively to the leading order in Nf/N c . We take 
3 cases for the background to help address the questions in stages, and gain intuition 
for how the solutions should behave. 

We first solve for the simplified system of a decompacitfied transverse coordinate 
of the D8 branes, which has been studied in its own right in [26], [27] (see figured] (a)). 

3 For this added complication, though, wc simplify the equations by using a perturbative ap- 
proach. In some sense, this is complimentary to smearing: one smears the branes to obtain non- 
linear ODEs to solve; we instead perturb the equations to obtain linear PDEs. 
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For this case we were able to find exact solutions of the partial differential equations. 
To our surprise we have found that whereas the solutions for the perturbations of 
some background fields behave, as we have expected, with a cusp at the location of 
the probe branes (a "A" shape), for other the behavior is of an inverted cusp dressed 
with a double hump structure (an "m", see figure [3]). We explore the ranges where 
we expect the supergravity to be a good description and find that u ^> ^/B? m , and 
that the perturbative results are good uptoM<^Cl/Q/ = Ani s /(g s Nf). Further, 
from this study, it becomes plausible that compactifying the £4 direction is possible, 
as the functions die off at large x 4 . 

Next, we address the compactified (but extremal) case both to view the effects 
of compactification, but also as a rough approximation of the "cigar" case at large 
u (see figure [H (b)). We first treat this case where we sum over the images of the 
uncomactified case, and then as a Fourier decomposition. From this we see that 
there is no issue with the dilaton tadpole constraint: the two cusps meets smoothly. 
We find that the decompactified limit emerges at large u. Both methods are applied 
because, although both series always converge, one series converges very quickly at 
large u, while the other converges quickly at small u. 

Finally we analyze the system of the near extremal D4 branes. We find that in 
this case the perturbation theory is good for u <C 1/Qf which is generically stronger 
than the u 3 -C R 3 D4 /g^ supergravity regime. Further, we find that the supergravity 
description is valid near u = Uk, with the additional constraint (Uk/ RdaY^Q f <^ 
l/£ s . This translates into the requirement that (T st / 'Mjj,) x / 2 A4-/V / / 'N c <C 1 where T st 
is the string tension, M g b is the typical glueball mass, and A4 = g\N c is the 't Hooft 
coupling. For large values of the radial direction the solution is obviously like that 
of the extremal compactified case. We use a Fourier decomposition to show that 
a finite expansion around the tip of the cigar is possible, and then implement this 
expansion for the first few terms. 

Once we have established the perturbative solutions, we proceed to analyze the 
stability of the system. We first show that the solution of the embedding of the flavor 
brane at the probe level persists also in the leading order backreaction. We further 
show that the fluctuations of the embedding, which correspond to scalar mesons 
in the dual gauge theory, are non tachyonic. Hence we shown that the system is 
stable at least for an action that is quadratic in the fluctuations. This is due to 
a cancelation between the electrostatic repulsion (CS action) and the gravitational 
attraction (DBI action). Hence, the above analogy with an electrostatic problem is 
not quite justified: the electric repulsion is canceled by a gravitational attraction. 
The only other force is that of the tension of the brane, which is restorative. The 
corrections to this force, while interesting, cannot change the qualitative feature of 
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Figure 1: The three cases studied: a) the uncompactified b) the compactified 

and c) the near extremal "cigar" case. 

stability while the perturbative analysis is valid (however, the effects of the non- 
perturbative backreaction is still an open question). 

The paper is organized as follows. In the next section (2) we briefly review the 
general setup of the problem, namely, the Sakai Sugimoto [10] model and the massive 
type II a supergravity action [28,29]. In section 3 we write the supergravity EOMs 
that incorporate the backreaction of the probe branes. We then introduce an ansatz 
for the metric which we substitute into the equations. The perturbative parameter is 
defined, and these equations are expanded. The gauge invariance, in the form of small 
coordinate transformations, of the system is discussed. In section 4 we present the 
solutions of the backreacted EOM. We start with the solutions for the uncompactified 
case, and then by summing over images the solutions for the compactified extremal 
case is constructed. We also use Fourier analysis to study this case. This enables 
us to determine the UV behavior of the near extremal case because the geometries 
are identical at large radius. The third step is to write down the solutions for the 
near extremal case in the region close to the horizon. In the following section we 
analyze the stability of the system. We first show that the solution of the EOMs that 
follow from the backreacted DBI+CS action are the same as those of the unperturbed 
solution. Finally we shown that the spectrum of fluctuations around this embedding 
is tachyon free and hence we conclude that to the leading order in Nf/N c the system 
is stable. 
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2 Brief review of the general setup 



Before we start the analysis of the backreaction of the flavor branes, we briefly 
review the two main ingredients of the general setup of the problem, namely, the 
Sakai Sugimoto model and the action of the massive type 11a supergravity. The 
reader familiar with these topics should skip to the next section. 



2.1 Sakai- Sugimoto (SS) model 

Constructing holographic models duals of gauge dynamics that admits confinement 
is by now a relatively easy task. Incorporating flavor chiral symmetry, on the other 
hand, turns out to be more complicated. A prototype model that includes both 
phenomena is the Sakai-Sugimoto model [10]. It is a model of a holographic dual 
for a 3 + 1 dimensional gauge theory with a continuous SU(Nf) x SU(Nf) flavor 
chiral symmetry which is spontaneously broken. It is based on the incorporation of 
Nf D8-branes and Nf anti-D8-branes into Witten's model [11]. The latter describes 
the near horizon limit of N c D4-branes, compactified on a circle of radius R (x^ = 
X4 + 2ttR x ) with anti-periodic boundary conditions for the fermions. The D8-branes 
are placed at x^ = and the anti-D8-branes at x^ = L. The gauge theory dual 
of this SUGRA setup is a 4 + 1 dimensional SU(N C ) maximally supersymmetric 
gauge theory, compactified on a circle with anti-periodic boundary conditions for 
the adjoint fermions, and coupled to Nf left-handed fermions in the fundamental 
representation of SU(N C ) localized at £4 = 0, and to Nf right-handed fermions in 
the fundamental representation localized at £4 = L. 

The basic assumption of the model is that in the limit of Nf <C N c one can ignore 
the back-reaction of the Nf D8 branes and Nf anti-D8 branes. As mentioned above 
the goal of the present work is to examine in details the back-reaction of the D8 and 
anti-D8 on the background. With the probe assumption the closed type IIA string 
background is given by : 



/ \ 3/2 / r> \ 3/2 

ds 2 = f ) \-dt 2 + Sydx'dx* + f(u)dxl] + ( — 
\ Rda J V u 



du 2 

+ u 2 dVt 2 , 



L/(« 



qp3 / „ \ 3/4 /TI \ 3 

^ = 6 -^Q„ e* = gj^-) , R 3 m = ng s N c ll /(«) = ! 



RD4 J V U 

(2.1) 

where t is the time direction and x l (i — 1, 2, 3) are the uncompactified world- volume 
coordinates of the D4 branes, X4 is a compactified direction of the D4-brane world- 
volume which is transverse to the probe D8 branes, dQ 2 is the metric of a unit 
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four-sphere and €4 is its volume form, and g s is related to the 4 + 1 dimensional 
gauge coupling by g\ = {2n) 2 g s l s . The sub manifold spanned by £4 and u has the 
topology of a cigar with u > Uk, and requiring that this has a non-singular geometry 
gives a relation between Uk and R x , 



2 (R 



3 \ 1/2 

L D4 



The parameters of this gauge theory, the five-dimensional gauge coupling g 5 , the 
low-energy four-dimensional gauge coupling g^, the glueball mass scale M g b, and the 
string tension T st are determined from the background (12. ip in the following form : 

2 / jj \ 1/2 -, 

t _L / 1 J_^V /2 J_^£ fo^ 

s * 27r/ 2V^-U=^ 27T/2 Vi? D4 ; ~ 27tt Rl 2WR%' 1 j 

where A5 = ^fiVc, M 9 6 is the typical scale of the glueball masses computed from 
the spectrum of excitations around (12. II) . and T st is the confining string tension in 
this model (given by the tension of a fundamental string stretched at u — uk where 
its energy is minimized). The gravity approximation is valid whenever A 5 ^> R x , 
otherwise the curvature at u ~ Uk becomes large. Note that as usual in gravity 
approximations of confining gauge theories, the string tension is much larger than 
the glueball mass scale in this limit. At very large values of u the dilaton becomes 
large, but this happens at values which are of order N^ 3 (in the large N c limit with 
fixed A5), so this will play no role in the large N c limit that we will be interested 
in. The Wilson line of this gauge theory (before putting in the D8-branes) admits 
an area law behavior [30], as can be easily seen using the conditions for confinement 
of [31]. 

The gauge theory dual to the SUGRA background (12.11) is in fact not four di- 
mensional even at energies lower than the Kaluza-Klein scale 1/R X since the masses 
of the glueballs are also M g j, = 1/R X , namely, there is no real separation between the 
confined four-dimensional fields and the higher Kaluza-Klein modes on the circle. 
As discussed in [11], in the opposite limit of A 5 R x , the theory approaches the 
3 + 1 dimensional pure Yang- Mills theory at energies small compared to 1/R X , since 
in this limit the scale of the mass gap is exponentially small compared to 1/R. 

The probe flavor D8-branes span the coordinates t,x l ,fl4, and trace a curve 
u{xi) in the (2:4, u) -plane. Near the boundary at u — > 00 we want to have Nf D8- 
branes localized at £4 = and Nf anti-D8-branes (or D8-branes with an opposite 
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orientation) localized at 24 = L. Naively one might think that the D8-branes and 
anti-D8-branes would go into the interior of the space and stay disconnected; however, 
these 8-branes do not have anywhere to end in the background (12. II) . so the form 
of u(x 4 ) must be such that the D8-branes smoothly connect to the anti-D8-branes 
(namely, u must go to infinity at x 4 = and at £4 = L, and du/dx^ must vanish 
at some minimal u coordinate u = uq). Such a configuration spontaneously breaks 
the chiral symmetry from the symmetry group which is visible at large u, U{Nj)l x 
U(Nf)n, to the diagonal U{Nf) symmetry. Thus, in this configuration the topology 
forces a breaking of the chiral symmetry. 

To determine the profile of flavor probe branes, one has to solve the equations of 
motion of that follow from the DBI + CS action that describes the probe branes. It 
is easy to check that the CS term in the D8-brane action does not affect the solution 
of the equations of motion. More precisely, the equation of motion of the gauge field 
has a classical solution of a vanishing gauge field, since the CS term includes terms 
of the form C 5 AFAF and C 3 AFAFAF. So, we are left only with the DBI action. 
The induced metric on the D8-branes is 

/ \ 3/2 / \ 3/2 

ds 2 DS = [-5-) [-dt 2 + ^jdx l dx j ] + (-£-} 



/(«) 



D \ 3 ,2 



u J f(u) 



R 



Di 



u 



3/2 

u 2 dQl (2.4) 



where u 1 = du/dx^. Substituting the determinant of the induced metric and the 
dilaton into the DBI action, we obtain (ignoring the factor of Nf which multiplies 
all the D8-brane actions that we will write) : 



S DBI = T 8 J dtd z xdx A d^e-*^- det(£) = y f dx 4 u 4 Jf(u) + (^^j yj^ 2 - 5 ) 

where g is the induced metric ( 12.41) and T§ includes the outcome of the integration 
over all the coordinates apart from dx±. The simplest way to solve the equation of 
motion is by using the Hamiltonian of the action (12.51) . which is conserved (indepen- 
dent of £4) : 

« 4 /M 



constant = Uq \/ f(u ), (2.6) 



where on the right-hand side of the equation we assumed that there is a point no 
where the profile ^(^4) of the brane has a minimum, u'(u = uq) = C@. We need to 



4 This type of analysis was done previously for Wilson line configurations. See, for instance. [30]. 
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(a) 



(b) 



Figure 2: The dominant configurations of the D8 and anti-D8 probe branes in the 
Sakai-Sugimoto model at zero temperature, which break the chiral symmetry. The 
same configurations will turn out to be relevant also at low temperatures. On the 
left a generic configuration with an asymptotic separation of L, that stretches down 
to a minimum at u = u , is drawn. The figure on the right describes the limiting 
antipodal case L = ttR x , where the branes connect at uq = Uk- 



find the solution in which as u goes to infinity, X4 goes to the values £4 = 0, L; this 
implies 

dxA = 2 J^ = L (2.7) 

with u' given (as a function of u) by (12.61) (note that u is a double-valued function 
of £4 in these configurations, leading to the factor of two in (12. 7p ). The form of this 
profile of the D8 brane is drawn in figure [2](a). Plugging in the value of u' from (12. 6p 
we find 

du , /2 r . 1 



dx A = 2 / — = 2R%1 I du- 



/(«(>)«o 

2 //? 3 \ 1/2 / r 1 z 1 / 2 

= P J l - y K / dz Z , (2.8) 

where i/k = uk/uq. Small values of L correspond to large values of uq. In this limit 
we have y K 1 leading to L oc \/R D4 /u . For general values of L the dependence 
of Mo on L is more complicated. 

There is a simple special case of the above solutions, which occurs when L = irR x , 
namely the D8-branes and anti-D8-branes lie at antipodal points of the circle. In this 
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case the solution for the branes is simply x±{u) = and 24 (u) = L = nR x , with the 
two branches meeting smoothly at the minimal value u = uq = Uk to join the D8- 
branes and the anti-D8-branes together. This type of antipodal solution is drawn 
in figure 12(b). It was shown in [10] that this classical configuration is stable, by 
analyzing small fluctuations around this configuration and checking that the energy 
density associated with them is non-negative. 

In general for L < nR x , there is a family of smooth configurations characterized 
by L or by the minimal value of u, uq. This class of solution is shown in [2(a) 

The Sakai-Sugimoto model has 3 dimensionful parameters : A5, L and R x , and 
gravity is reliable whenever A5 3> R x . The physics depends on the two dimensionless 
ratios of these two parameters; In the gravity limit the mass of the (low-spin) mesons 
is related to 1/L [13] while the mass of the (low-spin) glueballs is related to 1/R X . 
As discussed above, in the limit A 5 <C R x this theory approaches (large iV c ) QCD at 
low energies. This remains true also after adding the flavors, at least when L is of 
order R x . 

The thermal phases of the model where analyzed in [13, 14]. The back-reaction of 
the flavor brane at non zero temperature is not addressed in the present paper and 
will be described in a future publication. 

2.2 Massive type II A and 8 branes 

One expects p + 1 dimensional objects to naturally couple to a p + 1 form potential. 
Therefore, one expects a D8 brane to couple to a nine form potential. Conventional 
type IIA supergravity has no such form, and so some modification of the theory is 
necessary to describe the backreaction of D8 branes. This extension was first found 
by Romans [28] , and then further generalized to admit localized D8 solutions in [29] . 
The relevant kappa symmetric worldvolume actions were constructed in [32] . Further 
studies of D8 (-Dp) brane backgrounds (systems) are discussed in [33-38]. 

There exists another massive type IIA, constructed in [39]. This and the Romans' 
type IIA were shown to be the only "Higgs type" supersymmetric extensions of 
massless type IIA in [40], although a third was suggested in [39]. The massive type 
IIA given in [39] does not admit localized supersymmetric eight-branes, and so we 
restrict our attention to the theory of Romans [28,29] which we simply refer to as 
massive type IIA. 

The bosonic part of the action of the massive type IIA supergravity takes the 
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5m = / rf 10 xv^^f 2 + -i- / MF 10 

Zft 10 J Z Zfi -10 </ 

where • denotes contraction of indices with inverse metrics, and e is antisymmetric 
in indices and takes values ±1. In the above action F n+ i = dA n and 

F 2 =F 2 + MB 2 

F 4 = F i + ^MB 2 AB 2 (2.10) 
F 4 = F A -A 1 AH 3 + -MB 2 AB 2 . 

Li 

One notes that from the above definitions, F 2 may be absorbed completely by a shift 
in B 2 , but only when M/fl. One views this as a "Higgsing" where the degrees of 
freedom associated with F 2 become the longitudinal modes of a massive B 2 . The 
equation of motion for F 2 , therefore, must only be imposed in the massless limit. In 
appendix A, we include the equation of motion associated with Ai so that an M — > 
limit is clear. In the next section we turn to including sources in the action. 



3 Backreaction of D8 branes 

In this section, we will find the equations of motion that govern the D4-D8 systems 
of interest, including the contribution from the DBI + CS brane action. We present 
our ansatz, and the perturbative parameter we will use to linearize the equations, 
and finally present the separated linearized equations. Further, we find the small 
coordinate transformations that leave the form of our ansatz unchanged (to the 
order we are working): these are gauge transformations of the linearized equations. 

5 We use the notation of chapter 12 of Polchinski [41] 
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3.1 Finding the equations: ansatz and separation 

For the remainder of the paper, we will be concerned with D4-D8 systems, and 
because neither of these branes source (directly) either or B^, we set them to 
zero. After truncation to the = 0, B^ u = sector, the equations of motion for 
the type IIA massive supergravity are the following: 

e 2<P / 1 <A e 2 * 

R,u + 2V M V,0 = — [AF^^F ya2a . iaA - -g, v (F (4) ) J - -^g^M 2 

V Ql F aia2a3Q4 = (3.1) 
d^M = 
M = *F(i 0) 

U — fc - r ai«2"3 Q 4- r PlP2P3P4 ' 

Note that the equation of motion for A\ in the appendix is trivially satisfied. Again, 
one must only impose it's equation of motion in the massless limit. However, the 
equation of motion for B^ u imposes a constraint, arising from the Chern Simons 
term B A _F 4 A F 4 , on the four form (the last of the above equations). This constraint 
is easily satisfied for simple 4-form field strengths. Also, note that we have used 
the dilaton equation of motion (EOM) to eliminate R from the Einstein equation. 
This will be important below when we derive the equations when a brane source is 
present. 

We now turn to the modification of the equations of motion (13.1 ft by adding 

- K p (J d^e^^g~ p + / Vi) ( 3 - 2 ) 

to the action. Here we use g p to denote the pullback metric on the p+1 dimensional 
sub manifold defined by X a (£), 

dX"dX» 

and K is the appropriate constant involving the p-brane tension. In this action, we 
assume that it is consistent to set the world volume U(l) gauge field to zero, as well 
as ignoring any additional Chern Simons terms (which is appropriate for the cases 
we wish to consider). 
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There are two types of fields in this action: those that represent open string 
degrees of freedom (e.g. X Q (£)); and those representing closed string degrees of 
freedom (e.g. g a p). Of course when varying with respect to the closed string degrees 
of freedom, 10D delta functions appear, of which p+1 are integrated leaving behind 
a (10 — (p + 1)) dimensional delta function source term, as we should expect. Varying 
the above action with respect to the p + 1 form potential adds a delta function source 
to the form fields equation of motion of the generic form 

Note that the product e w ■ e p+ i is sensitive to the orientation of the submanifold de- 
fined by X a . For example, in the case of the SS model with the antipodal embedding 
of the D8 branes, there is both a positive delta function and a negative delta function 
in accounting for the orientation reversal of the brane (it is oriented in the +u 
direction). For D8 branes, however, one should replace d* F with dM because M is 
the term that appears with F\q in the action. Hence, M is in fact piecewise constant 
in backgrounds with localized D8, as we will see below. 

We restrict ourselves to embedding functions of the form X a {^) = £ a and the 
remaining X 1 are arbitrary constants. Varying the full action with respect to the 
dilaton and graviton is now straightforward, and the equations of motion are 

R, v + 2V M V„0 - — (4F M a ^F ua2a3a4 - -g, v (F {4) ) J + — 9flu M 2 
Ks24o , >Ffel<w f^dX-dX" 1 \ _ 

+ g 6 ~7=| I 9 P Q^a Q£b 9an9fi» ~ 2^ J 6 i X ~ X \ X >) ~ U 

r + v%v y - 4<r (v,0) (v,0) - ^k e ^ y^jk^L s (x a - X°{x)) = 

-dM - K 8 6 (x a - X a (x)) (£l ° ' ' £9) = 



(2«? ) ' W7 9! 

*F {10) = M. (3.5) 

with other equations of motion left unchanged. The delta functions appearing above 
may be simplified by taking them to be functions only of £4, 5 (x a — X a (x)) = A(x^), 
which is appropriate for the antipodal embedding in the Sakai Sugimoto model. As 
expected, only the RR couplings to the branes are sensitive to the orientation of the 
branes. Further, the epsilons appearing above take values ±1, and do not contain 
factors of J—g. lI 



3 As a simple check of the above signs and numerical factors, one can simply check the follow- 



11 



The tensor structure of the Einstein equations can be easily read: the delta 
function strength is proportional to the metric and dilaton, and comes with a + sign 
for a direction along the D8 brane, and comes with a — sign for those not along 
the D8 brane. Although above we have written the effect of a D8 brane, the above 
arguments work also for an arbitrary p brane: it simply changes which RR form field 
equation of motion gets the delta function source, and how many directions of the 
Einstein's equations get (— ) vs. (+) delta functions. 

For the remainder of this work, we will take the solution to the M and Fiq 
equations of motion to be 

M = ±2.^/2 = ^^ 
*F 10 = M (3.6) 

where the + is used on one side of the D8, and — is used on the other. 

The remainder of the paper will be devoted to solving the remaining equations 
of motion perturbatively, and the perturbative control parameter will be explained 
shortly. However, at this point an important note is in order: delta functions in 
codimension 10 — (p + 1) ^ 1 have singularities at the location of the delta function. 
Codimension one is special in that the Green's function is of the form \x\ 0, and 
hence is finite at the source. We therefore expect that the perturbative approach is 
most natural for D8 branes, as the back reaction can be made small, even close to 
the brane. Hence, the perturbative approach that we take may not be as natural for 
higher codimension branes. 

To define the small parameter in our expansion, we make the following obser- 
vations. Given the solution of the M, F\q sector, the new (relative to the mass- 
less IIA equations) terms in the equations of motion come with the powers of 



ing. The coefficient in front of g C p b d a X a dbX 13 g all gp v can be checked against that of R^ v . Before 
using the dilaton equation of motion, the Einstein equations contain ~^Rg^ Ul and this coeffi- 
cient must match that of g'£'d a X a dbXP g aiL gp Vl except that one multiplies the latter by an addi- 
tional — A"g2K^Q-y/— g v f \J—g. This is because they are obtained from similar terms in the action, 

respectively. The factor in the dilaton equation is also obtained 

similarly, as the coefficient of R is obtained by varying -^-J— ge~ 2 ^R w.r.t. <fi and the delta func- 

tion coefficient is obtained by varying —K sy /—ge~^, and so a factor of 1/2 arises. The second part 
of the delta function term in Einstein's equations is similarly found by checking that one adds \g^v 
times that of the dilaton term. 

7 This terminology only makes sense for branes of dimension D — 1, as such branes split the space 
into disjoint regions. 

8 we refer to this behavior as a cusp 
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2k\qK^ ~ g s Nf. This is what we shall use as a control parameter for our per- 
turbative expansion. From now on, we will simply take 2k\ q K%/2 = Qf/g s as the 
definition of our small parameter Qf. Another way of phrasing this is that in the holo- 
graphic limit, there is a scale R c such that g s ~ (R c /l s )/N c where R c /l s ^> 1 is large 
but held fixed. Hence, one may view our perturbation as limit on g s Nf ~ Nf/N c , 
which is the basis for the probe approximation. 

To solve the equations, we will still need to take an ansatz, and we motivate it 
as follows. The eight brane doesn't directly couple to F A , and further the 5*0(5) 
symmetry of the 4 sphere is not broken for this brane configuration. Hence, we 
take that there is no change in F A to leading order in Qf. Further, because we take 
the solution M oc ±Qf/g s , so that in the Einstein equations, the term g^ u M 2 e 2 ^ ~ 
O (Q}) may be ignored. 

Therefore, we assume that only metric and dilaton perturbations are necessary, 
so we take a general ansatz of the form 

ds 2 = e 2A(u, Xi ) (_ df 2 + dx T^ + e 2B(u, X4 ) dx 2 + ^2 + e 2C M dQ 2 

(f)(u, x A ) = x A ) + 2A(u, x A ) + 2C(u, x A ) (3.7) 

where Q c = 3R\/ g s = SnNJl, and is the volume form of the unit four sphere. It 
is clear that the F A equations are trivially satisfied: dF A = because Q A is closed, 
and d * F A = because *F A is some function of £4 and u times du A dx A Aj dxi, and 
is therefore closed. 

In the above, we will expand the above functions as 

A(u,x A ) = A (u) + QfA^u^i) B(u,x A ) = B (u) + QfBi(u,x A ) 
C(u 1 x A )=C (u) + Q J Ci(u ) x A ) G(u ) xA)=G () (u) + Q f Gi(u ) x A ) (3.8) 

4>{u, x A ) = 4> (u) + Qf(j)l(u, X A ) 

where the subscripted functions are solutions of the Qf = equations. We linearize 
and explain how to separate them for the Sakai Sugimoto model in appendix B, and 
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summarize the results here. One must solve the decoupled system 



u(u 3 -U 3 K ) (u 3 -U 3 K ) 2 
54uF x , 2u{Q c g s )lVZ 



(u 3 



U k) ^-Ul) \J(Q c9s )*(l- U Jt) 



A = (3.9) 



3^F 2 + 



3(Au 3 - U 3 K )d u F 2 < g a Q c u 3 ^F 2 



u(u 3 — Uk) (u : 
4u(Q c g s )3 a/3 / w 



it- 



A = 



(3.10) 



_ 49 2 7 4u 3 Q c <? s ^ 4 0i 2{u 3 -WDdufa 36^0! 



i 



3( U 3-^) 2 M^-t/l) (u 3 -U 3 K y 

i>(Qcg s )%V3 



+ / 7 ry A = (3-11) 



where 



I 5 (x 4 ) + 5 (x 4 — irR x ) if x 4 = x 4 + 2^7?^ 
] 5 (x 4 ) if x 4 non-compact 

and then identify the physical degrees of freedom 



(3.12) 



3 - 

A x = — Fx + — F 2 0! 

5 10 1(T 

Cj = — Fi + -F 2 - —0i 
1 10 5 10 r 

1 1 1 - 

B 1 = G 1 = --Fi{u,x A ) + -F 2 (u,x 4 ) - — 0i (14,0:4) 
5 5 10 

--^0^X4) + - (u ,_ E/ 3 ) (3.13) 
0i = -0i + 2Ai + 2Ci. 
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Above, we have made the obvious notation that <pi is the first order correction to 
the physical dilaton. One may read off the combined solution by plugging in these 
to the equations (IB. 71) in the appendix. 



B x ^ B L 



3.2 Gauge Freedom 

Here we identify the gauge (coordinate transformation) freedom as those transforma- 
tions in u and x± that leave the metric diagonal (preserves the form of our ansatz). 
Indeed, 

2 2 5A(h,x 4 ) 

2 u 

3 X(u, X4) 

4 u 

3X(u,X4) 3 U^X(u, Xi) 
4 u + 2u(u 3 - U%) 
1 f Qc9s(d X4 d X4 X(u,x A )) 

"3 / 7 — ^ — ( } 

3 AM 3 U z K X(u,x 4 ) 

1 A(M,g 4 ) 

C-i — > Oi + 

4 it 

leaves all equations of motion unchanged, and is exactly a coordinate change in u 
and X4, namely 

u — » it + A(it, £4) (3.15) 
1 f Q c g s {d XA X (u,x A )) 

u 3 (l-- 



£4 — > X4 / « U 

3./ o I/^ 2 



Hence, one of the degrees of freedom above is pure gauge. However, there is an added 
complication. If we eliminates <pi using such a gauge transformation, the cusp in <fri 
generates a delta function in the gauge transformation for Bi, and hence Bi is no 
longer a smooth function: it contains a delta function. Hence, we conclude that for 
the unsourced equations one may choose which degree of freedom to eliminate, but in 
the sourced equations, only B\ may be eliminated. However, as shown in appendix 
B, it is more convenient to not eliminate B\ completely, but rather to take B\ = G± 
as the choice. 
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Given the above transformations, we can immediately see that F\ and F2 of the 
last subsection are gauge independent. The remaining gauge dependent quantities 
0i, B\ and G\ do not admit a gauge independent combination. Further, given the 
equations (I3.13p . only the equation B\ = G\ is not gauge covariant. Therefore, we 
will sometimes refer to this as a gauge fixing. 

4 Solutions: the linearized backreaction 

Here we will analyze the differential equations of the last section in three separate 
cases: 

1. Uk = decompactification limit: In this case we take Uk = 0, in a limiting 
sense of the background. In this limit R\ oc R^ d ^/Uk becomes infinite, and so 
£4 decompactifies. 

2. Uk = 0, 24 compactified: In this case, we note that while the Uk = limit 
has decompactified X4, one still has the isometry 24 — > X4+ constant. Hence, 
one may orbifold by this isometry and compactify £4. We will parameterize 
this compactification using the same radius, R 2 X = (A/^R^^/Uk- One way 
to think of this parametrization is that we have taken the spacetime to be 
that of Uk = while still requiring that x 4 is compactified: we choose to 
parameterize the compactification of £4 such that we may compare easily to 
the Uk 7^ case. In this way, we have taken the spacetime to be the cylinder 
to which the cigar asymptotes, and so this analysis gives the u ^> Uk behavior 
of the Uk 7^ case. In all compact cases, we will be considering the 
antipodal embedding, L = 7iR x , which for concreteness we parameterize by the 
embedding a; 4 = 0, nR x . 

3. Uk 7^ 0: In this case, we analyze the equations as is. We make some basic 
comments about the nature of the Fourier transformed equations, and note 
that the point u = Uk is a regular singular point, and hence a finite convergent 
series about this point exists. We expand the solution about the tip of the 
cigar. 
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4.1 Uk — decompactification limit 

In the Uk — > limit, the differential equations become 



3 ^ 1+ H a „ Fl _^ Fl+ % alFl+ w^ = o 

3^ + -^ 2 + %^F 2 - 4 ( 3Q f)' A = (4.1) 



u w 



2 ~ 2 . 4Q C <^ 2 ~ 4(3Q c g s )2A 

- - ^-r^ + 3m§ = 0. 

In the above equations, we take all functions to be functions of the form 



Fi(u,x 4 ) = uKi(q) q 

{Qcg s ) 2 

4>i(u,xA) = uK 3 (q). (4.2) 

This has the effect of changing the delta function in £4 into a delta function in q as 
A (#4) = A'(q)y/u] [QcQs]- Further, we take just a single brane so that A' = 5(q). 
Taking the resulting equations, and multiplying them by u, we obtain ODE's with 
delta function sources 

{^q 2 + d 2 q K x + ^qd q K x - A2K X + 2V35(q) = 
Qg 2 + d 2 q K 2 + ™qd q K 2 + 12K 2 - 4V35(q) = (4.3) 
(y 2 + l) d 2 q K 3 + 3qd g K 3 + ^K 3 - V35(q) = 
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One constructs the delta function solution from the vacuum solution. The vacuum 
solutions to these equations may be written 

1 1 2 . .4 



(3402g 13 + 22113g n + 57915g 9 + 77220g 7 + 54054g 5 + 18018g 3 + 2002g) 
+K b 



Ko = K, 



(3g 2 + 4)i 

2 -fq In ( v/V + 4 + gV3) + VW+IVS (g 6 + f g 4 + f g 2 - if) 

(3g 2 + 4)i 

g In ( - qVs) ~ VW+^V3 (g 6 + f g 4 + f g 2 - if) 



(3g 2 + 4)2 



K * = K *WTl + K iWTl (44) 

To obtain even (in g — > — g) convergent (for large g) quantities with cusps, we may 
construct the combinations 



-256^/ 

(3402g 13 + 22113g n + 57915g 9 + 77220g 7 + 54054g 5 + 18018g 3 + 2002g) 



3402 / 1 1 2 4 

* = ^ooi-^ir(i2- + / +? 



(3g 2 + 4)f 



K 2 = 2 10 v / 3 ^-^ (4.5) 

f3o 2 + 4)l 



(3g 2 + 4)1 

^ ln (^S^vl) + 2 V/^TIVS (g 6 + f g 4 + f g 2 - if) 



2 7T~n ITs 



(3g 2 + 4) 



2 



K 3 = 2V3^^ + N 3: 



3g 2 + 4 d 3g 2 + 4 



The above have been written with the cusp solution first, and then an even function 
that converges (with coefficients iVj). We have not been able to determine physical 
boundary conditions that fix iVj, and so we will leave them arbitrary when possible. 

To graph them, however, we take Ni = and show these in figure [3j Further, 
note that the function K 3 has a larger characteristic width, as it only converges as 
1/g. This will be important when we compactify x^. 
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(a) K^q) 



(b) K 2 {q) 




(c) K 3 (q) 

Figure 3: Graphs of Ki(q), setting N 2 = N 3 = 0. 

The height of the above functions grows as u because the peak happens at a fixed 
value of q, giving just a constant contribution times the dressing factor of u. Thus, 
one expects the perturbative approach to be valid up to u <C l/Qf- This will be 
generic for later sections as well, as the decompactified behavior emerges at large u 
in the following sections. 

We also wish to characterize the width of the "spike" in each graph. One way 
is to make sure that variations happen on scales larger than string scale. The slope 
of the graphs is largest in the vicinity of the spike, and this slope is determined 
by its q = behavior, and is therefore just a constant. Therefore, the x\ slope is 
simply d^i/dx^ = QfU 3 / 2 R~^/ 2 x constant. The physical length that this corresponds 
to, however, is ds = (u / Rd^I 4 ^ dx±, and we require that dfa/ds -C l/£ s . This gives 
u 3 < RdJ(Q}^) RhJigtNf). Recalling the conditions above u < l/Qf and 
N c ^> Nf, this condition follows, and so is not a new piece of information. 

One may also wish to characterize the width when the linear part is no longer the 
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dominant, and so characterize the width of when other "features" become important. 
This occurs when the q coordinate becomes order 1, and so translates into X4 oc 
R^l/u 1 / 2 . Again, translating this into a physical distance we find s = R 3 ^u 1 ^ 4 ^> £ s , 
where we have required that this distance be greater than string scale. This gives 
a lower bound on u, however, it is the same lower bound coming from the Ricci 
scalar R oc l/fw-R 3 ) 1 / 2 C l/fj for the supergravity approximation. We see that we 
trust the supergravity to describe the backreaction above u 3> l A s /R? Di) and that the 
perturbative results are good up to u <C 

Of course one may take the last two constraints on u and turn them into a unitless 
constraint on the parameters. We find that this is l/Qf l\jR? D/ ± — > g 2 Nf/N c <C 1 
which we can see is weaker than the other constraints Nf/N c < 1, g s < 1. 



4.2 Uk = with £4 compact ified 

Here we will examine the Uk = case with x 4 compactified as explained at the 
beginning of section 4. However, a few brief words are in order. We will do this case in 
two ways: by summing the images from the last section, and by Fourier decomposing 
them. These two approaches are complimentary because one expects the sum on 
images to converge quickly for large u (when the images are well separated), and as 
we will see, the Fourier modes converge extremely quickly as u — > 0. Further, we 
will be considering the antipodal embedding, L = 7rR x , which for concreteness we 
parameterize by the embedding x 4 = 0, nR x . 



Sum images of Decompacitification 

Looking at the Fourier transform (below), there is no problem with compactifying 
a?4, however the solution for K 3 above seems to have a problem. Note that if we 
take the above K 3 and sum over images in q for some fixed value of u, we expect a 
divergence, as the function only converges as 1/q. However, switching back to the 
u, £4 language, we find that the n th image of K% is 

3 

k,t = ra > ~ V " V | " (4-6) 

*(& + ™) + 9 ^ ' 

where the right hand side is it the large n behavior. This, however, is actually a 
homogeneous solution to the original differential equation we started with. This 



X4 



irn 
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suggests a solution to this difficulty, and we take instead 



3^ 



+ Tin 

tlx 



9 



n\n + C n \ 



(4.7) 



where this is just adding a zero mode of the differential equation. The sum of this 
function in n converges, as the large n behavior is order 1/n 2 . 

Now we make one more final comment. To fix C n , we require that 
lim^oo 0! (w, x 4 = ttR x /2) = 0. This is roughly requiring that when you are as far 
away from the branes as possible that the perturbation should be 0. The solution to 
this constraint is C n = 1/2. Of course a different set of C n could be chosen in such 
a way as to not affect the sum: this, by definition, is unphysical, as none of the field 
values would change. 

Therefore, for the compactified case, we take 



/ 



00 

E 

n=— oo 



3 u 



■'4 



Tin 



\ 



U K 



tlx 



9 



7r n + 



(4.8) 



The other functions we take N 2 = and simply sum on images. We plot these in 
figure HI 

In figure H] we have summed many fewer modes for Fi because this function 
converges very fast (as 1/g 14 ), and so fewer images are needed. Further, we can see 
that one must go to much larger values of u to obtain well separated features for 
One could have guessed this from the results of the decompactified case: the 



a- 



characteristic width of the functions Ks(q) (associated with 0i) is much larger than 
for the other functions K\ and K 2 . However, one can see the shape of the functions 
approaching the decompactified case: F-y has a single cusp like —\x^\ and 0i and F 2 
have cusps like \x±\ surrounded by two maxima. We show plots for various u in the 
next subsection where we consider the Fourier transform of the equations. 
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(a) F^x^/Uf 



(b) F 2 (u,x 4 )/U 



K 




(c) (j>i(u,X4,)/U K 

Figure 4: F\/Uk as a function of x^/R x for u/Uk = 100 using images from n = 
—50 ■••50, F 2 /Uk as a function oixi/R x for u/Uk = 100, n = —500 •••500 and 
4>i/Uk as a function of x^/R x for u/Uk = 10, 000, n = —500 • • • 500. 



Fourier decomposition 

We start with the same equations as the last subsection 



3^i 



3dlF 2 



128^ , gsQcdl^ hAF 1 2{Q c g s )l^> 



+ 



+ 



u 



if 



u- 



12d u F 2 g s Q c d 2 XA F 2 A(Q c g s )IV3 



if 



li- 



lt 



a. 



(Q c g s 



(Q c g s ) 
-A = o 



-A = 



4Q c g s d 2 X4 t 



2d u fa , HQcgs)^Vs 



if 



u 



If 



II 



(Q c g s 



—A = 



(4.9) 
(4.10) 
(4.11) 
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and Fourier transform them 



39lFm , + IgMki _ _ + ih^S +,(-') i = o ( 4.i2) 

i ' i/O It — rt — 



38 ; Fm2 + _ _ "^(1 +, = Q (4.13) 



-4^0 m>1 + " + ^^^^ = o. (4.14) 

where we have assumed that £4 is periodic with 

x 4 = x 4 + 2vr^, = __ = -_ (4.15) 

as in the last section (the reason for this parameterization is explained at the begin- 
ning of section 4). 

The above equations can be brought to a simple and familiar form by the following 
change of coordinates and functions 



GmA 



FmAu) 



F m o(u) 



9m 2 U K 



3 

M2 



G 



m,2 



m,2W - r 

U2 



i(«) = G m3 I \l 9m ^ I U2 



it 



/; 2 



The above differential equations become 

p 2 d 2 p G m>l + p«9 p GW + (V - 9 2 ) G m>1 = _ 972(l + (-ir)m 4 ^ (41g) 



2?rp 

H 
2vrp 

1 

2tt~p 



p 2 d 2 p G rn , 2 + p0 p G m , 2 + ( V - 3 2 ) G mil = 1944(1 + ( (4.17) 
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These then are of two forms: modified Bessel equations with a simple monomial 
soured, and an exponential function with a simple monomial source. The general 
solution to the above equations is [42] 

G m>1 = CJ 9 {p) + C 2 K 9 (p) 27rp4( _ 13)(5) lF2 ^-T'2'Tj (419) 

C C T (n) I CK(n\ I 1944(1 + (-1)^4 / 5 iy\ 

G m3 = C,e* + C 6 e^ + 3(1 + Linm P) ± ^ ^ ZMM (4.21) 

where / and K are the modified Bessel functions, \F 2 is the generalized hypergeo- 
metric function 0, and Ei(a, x) is the exponential integral function. 

/■oo 

Ei(a,x) = / e { - yx) y- a dy (4.22) 



and we have used Re (Ei(l, —p)) = Ei(l, —p) +7ri to remove the — ni associated with 
going around the branch point at x — 0. 

Here, one may worry about the convergence of the above Fourier decomposition 
because the coefficient above depend on m in positive powers. However, recall that 
we wish to sum on m for fixed u, not fixed p. In fact the factor of m completely 
cancels out of the above coefficients once returning to the u coordinate (p 4 oc m A /u 2 ). 
The only m dependence comes about in the arguments of the homogeneous and non 
homogeneous terms. Therefore, we may effectively analyze convergence of the Fourier 
modes as convergence in the variable p — >• oo, as this is the limit to which m — ► oo 
corresponds. The inhomogeneous solution to (I4.2ip is indeed convergent and admits 
a power series expansion about p = oo. Therefore one does not wish to turn on the 
growing exponential (this would not converge summing on m), and the shrinking 
exponential is simply negligible. Hence, we may set C5 = C$ = for all m and get a 
convergent series. 

The remaining equations, however, deserve some special treatment. The bessel 
equation (and the equation for 1F2) have an essential singularity at p = 00. There- 



solutions to these equations are known as (modified) Lommel functions. However they are 
related to generalized hypcrgcomctric functions. We opt to use the notation of hypergeometric 
series, as these are more general, and perhaps more familiar to the reader. 
10 we do in fact mean 1F2, not 2F1 
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fore we consider the asymptotics of the above functions for large p: 

Up) - -k=f 



and [43] 



L F 2 (a 1 ;6 1 ,6 2 ;x) -> (4.24) 



r(6i)r(6 2 ) , ^i (ai _ fcl _ &2) 



-x)^ ai " t>1 ~ t>2 ^cos - Ui-foi-6 



7T 



0rr(ai) v ~ y ~" V 2 V~* " X 2 ' 2 

After substituting in x = p 2 /2 and choosing appropriate branches for the square 
roots ((— l) 1 / 2 = one finds 



P 4 2772p5 



i^lzllMj , v^^ ^ (425) 
P 4 60p^ 



We therefore use the following combinations 

Gm , = C ^ 9( „) _ 972(1 + (-in^| / ^jjig _ 

9W 2vr(-13)(5) 1 p 4 1386 I V ; 

Gm , = CAW + 1944(1 + (-1)")^ ( 4jM) + ^ (4 . 27) 
W 2tt(-7)(-1) \ p 4 30 I V ; 

G m3 = C«> + 3(1 + '- 1 > m >'^ Ei ' 1 - g> + t^' 1 ' =glM) . ( 4. 28 ) 

27T 

It is now a simple matter to replace the definition of p above and sum the series. 
Although we offer no analytic proof here, the above functions can be seen to converge 
quickly enough for large m. We note that because the above functions are functions 
of m 2 /u, the convergence in m and u — > are connected. As promised, the Fourier 
expansion converges more quickly for smaller u. 

In the last two equations, it should be noted that for large u, the particular 
solution selected dominates over the remaining inhomogeneous solution. However, 
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in the first equation, the homogeneous solution dominates. In the following, we still 
set C 2 = C 4 = C 6 = 0. 

Solving the m = case is trivial, but for completeness, we give the solutions for 
this as well 

-^0,1 — ZFZ <~ ^0,5^ + — g" 



; 0,1 



657T u 

24ukll Co 

7tt u 

1 i i i 

3£/j?- In(w)u2 Qu^Uk 



F , 2 = ^ + -X + Co,4 (4.29) 



7T 7T 



+ 2C 0A u2 + C ,2 (4.30) 



and we again set the unfixed constants above to 0. 

We show here the plots of the physical fields (f>i,A\,Bi = G±, C\ for the first 150 
modes in figure [51 
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(c) B l (u, x 4 )/U K (d) d {u, x 4 )/U K 



Figure 5: Graphs of <Pi/Uk, Ai/Uk, B\/Uk = Gi/Uk, C\/Uk graphed as a function 
of X4/R x . Each graph represents a different u value of u = 1, 5, 10, 20, 40 in ROYGB 
order. All graphs use the first 150 modes (m = 0...150), of which only even m are 
non zero. 

We may wish to ask what the long distance behavior is in u/Uk- For this, we 
plot the long distance behavior of the independent modes Fi, F 2 , 4>\ in figure [61 
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(a) F^x^/Un 



(b) F 2 (u,x 4 )/U 



K 




(c) 4>i(u,X4)/U K 

Figure 6: F\/Uk as a function of x±/R x for u/Uk = 100, F 2 /Uk as a function 
of 24/ R x for u/t/ft: = 100 and (pi/Ux as a function of x^/R x for u/Uk = 10,000 
(bottom) all graphed using the first 400 modes (m = 0...400). 



The graphs in figure [6] are directly comparable to the plots in figure IH however 
there may be discrepancies between the two. For example, in figure one notes that 
the function does not quite touch at x 4 = nR x /2 as it is required to for figure 
HI This is merely a matter of setting the correct zero mode for the homogenous 
equation in this case. However, to make the plots match exactly in all cases, one 
may, if one wishes, Poisson resum the results from the last subsection, which will fix 
all constants above in terms of the constants from the last subsection. We, however, 
have not done this. 

We should note also that the above analysis allows for more generality that sum- 
ming the images. In the decompactification case a very specific functional form was 
taken, to which many homogenous solutions do not conform. This property is in- 
herited when summing on images. By the orthogonality and completeness of the 
trigonometric functions, we are guaranteed to generate all solutions to the original 
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gravitational ansatz when using the Fourier analysis. 



4.3 U K ^ 

Fourier decomposition 

Fourier decomposing the general equations at the beginning of section 3 gives 

3(4 M 3 - U 3 K )d u F m>1 m 2 ^u 3 F rnA 



u(u 3 - U K ) 



[u 3 -U K ) 2ti(u 




+ 9Mv Sv ra ;' ' ./^w = o (4.31; 



W 2 p i 3(4« 3 - U 3 K )d,F m , 2 m^n'F 2 
jC,i m .2 H 



w(w 3 - U K ) 




-Ad: 



lguyg + (-l) m ) 
2vr(u 3 - 



2 (n 3 - 7U%) dj m>1 _ 36uUl4> m>1 4 m 2 « 3 ^i, 1 
^ ~ _ « (^ 3 - Ufc) (u 3 - t/i) 2 3 („3 _ [7^)2 
6uv^(l + (~l) m ) / u 



2tt(m 3 - U 3 K 



where now the periodicity of 




(4.32) 



0. (4.33) 



x 4 = X4 + 27r^, = __ = -_ (4.34) 

is dictated by smoothness of the gravitational solution near the tip of the cigar 
(u —* Uk)- Recall that we will be considering the antipodal embedding L = irR x , so 
that the delta functions are located at £4 = 0, ttR x . 

It is clear from the above that these functions are functions only of u/Uk as Uk 
is the only dimensionful parameter left. Under the coordinate change u — ^- and 
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the redefinition of fields Fi = Fj x Uk and <pi = (pi x [/# we find 

3(4u 3 -l)^ mil 27m 2 u 3 F mil 



u(u 3 - 1) 4 (u 3 - l) 2 

^ + 9a»(i + (-ir) =0 

(u 3 -l) 2vt(m 3 -1)I 



3<9 2 F m>2 + 



3(4£ 3 - l)d a F ma 27 m 2 u 3 F m , 2 
4 (w 3 -l) 2 



u{u 3 — 1) 
18u 3 (l + (-l) m ) 
2tt(u 3 - 1)1 



(4.36) 



2 ~ 2 (m 3 - 7) du4> m ,i 36u0 mi i n m 2 M 3 m ,i 

-4^0 m ,i — 3 — -2 + 9— -2 

u(u 6 -l) (u 3 - 1) (u 3 - 1) 

6u 3 (l + (-l) m ) , „, 

+ i 1 >> = o. 4.37 

2tt(m 3 - 1)1 

The above differential equations are difficult to solve, even excluding the inhomo- 
geneous piece. The homogeneous parts of the equations can be seen to have 5 regular 
singular points at u — 0, u; , u 1 , lu 2 , oo where uj = — | + &i is a third root of unity. 
However, as all of the singularities are regular, one may go about finding a Laurent 
series expansion about any given point (using standard textbook techniques). Rather 
than doing this in the Fourier basis, we find that it is easier to leave the functions of 
£4 intact, and expand the function about a specific u value, and then solve for the 
functions of £4 multiplying each power individually. 

u — > Uk expansion 

First, for expanding around the point u = Uk, and to make contact with much of 
the literature, we change to the coordinate z defined by 

u 3 = U\ + U K z 2 . (4.38) 
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In these coordinates, the (non algebraic) equations read 

QMUl + z^F, 27 {U 2 K + z 2 )d 2 z F 1 27 (U 2 K + ?>z 2 )d z F 1 
{Ujc + U K z 2 )h 4 U K 4" U K z 2 ' 4~ U K z 3 



54F X 2(3Q c g s )U(U K (U 2 K + z 2 )^A 

+ -U^ = ° ( 430 > 



Q c g s {U 2 K + ^)d 2 Xi F 2 27 {U* + z 2 )d 2 z F 2 27 {U% + 3z 2 )d z F 2 
{Ujc + U K z 2 )h 4 U K 4" U K z 2 ' 4" U K z 3 



4(3Q c g s )^(U K (U 2 K + z 2 )^A 

= (4.40) 



U K z 3 



4 Q c9s (U 2 K + ^)d 2 X4 h _ 9 (Uj c + z 2 )d 2 J 1 _ 2 {2z 2 -W 2 K )d z ^ 
3 {Ujc + U K z 2 )h 4 U K 1 U K z 2 U K z 3 



z A U K z 3 



B 1 =G 1 = --F, + -F 2 - 3 d Mu 2 K + * 2 ) + 1 MQU 2 K - z 2 ) 



5 5 5 z 10 z 



2 



(4.42) 



after a bit of cleaning up. 

Now, to expand about z = 0, we expand the functions multiplying each differen- 
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tial operator that acts on the fields F{, <p above. 



Qc9sd 2 x ,F x + 27 Uk%F\ + 27 U K d z F x _ 5AF\ + 2(3Q c g 8 )?y/U^A = Q 

z 4 A z 2 A z 3 U K z 2 z 3 

Q c g s d 2 Xi F 2 | 21 U K d 2 z F 2 | 27 U K d z F 2 A(3Q c g s ^VU^A _ Q 

z 4 4 z 2 4 z 3 z 3 

AQ c g s d 2 Xi fa 9U K d 2 Jt 27U 2 K d z fa 3GU K fa , f (3Q C *? S )5 ^A 



z 2 z 3 z 4 z 3 



(4.45) 



B 1 = G 1 = --Fx + -F 2 - - x zyi + --^i 4.46 

5 5 5 z 5 z 2 



We will now assume that all of the Fourier modes of the above fields have expansions 
about z = with a finite number of negative powers. This allows us to simply count 
powers in z and neglect any terms that are not of leading order. For example, this 
allows us to drop the term F\j z 2 above. The rest of the equation is homogeneous in 
z only if F\ ~ z as z — ► 0. Likewise, we conclude that to leading order Ft = zLi(x±) 
and 0i = zL 3 (x4). To leading order in z, the above equations become 



27/7 / 3/7 

d x ,d x M^) + -r^L^xA) + 2<-^A = (4.47) 
4Qc9s V Qc9s 



07TT / Q/7" 

^ 4 <9 X4 L 2 (x4) + -r^L 2 (xA) - 4W-^A = (4.48) 



07TT / qtt 

<9 X4 «9 X4 L 3 (*4) + -r^L^xA) - J — ^A = (4.49) 

B 1 = G 1 = z\ (L 2 - L0 . (4.50) 
5 
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These are easily solved. 



Li(x 4 ) 
L 2 (x 4 ) 
L 3 (x 4 ) 
R 2 , 



+ 3 
2 

+ 3 



&Qc9s 



sm 



sin 



sm 



£4 

£4 

Rx 



27U, 



K 



(4.51) 
(4.52) 
(4.53) 
(4.54) 



where R x is the periodicity of x 4 = x 4 + 27ri? x defined before. Note that this is 
exactly the kind of behavior one would expect because z\sui{xa/ R x )\ = \y±\ when 
switching to "cartesian" y 2 = z cos (x±/R x ), y\ = z sin {xa/ R x ) coordinates, and is 
similar qualitatively to a D8 in flat space. 

One may continue this process and in fact get the above functions to the next 
order. From the general differential equation in z, one may easily read that the 
expansion will be in terms of odd powers of z. This is because the functions being 
expanded are all {U\ + z 2 ) n , and so coefficients of even (odd) powers only mix with 
coefficients of other even (odd) powers. The next term, therefore, should be of order 
z 3 . We solve the resulting equations, and find 



Fi (2,2:4) 
F 2 (z,x 4 ) 



£4 



2 

— z sin 
3 \ R x 



17 



^(3sin( W -sin/'^M 



R, 



R, 



2 x -zsin 

3 V R x 



108 
14 



I J 2 



sm 



3|x 4 | 
Rx 



108 



U 2 



1 zs in( l -^)+0(z b ). 
Rx 



+ 0(z 5 ) 

(4.55) 



One may analyze the above functions for the various length scales of these cusps. 
One finds that they slope is directly limited by dfa/ds = d<pi/ dy\{UK / Rda)^^ ~ 
(U K /R D4 ) 3/4 Q.f < 1/4. This gives {ggT^Nf « (T st /M 2 gb ) l / 2 \ 4 N f /N c < 1. 

Further, the onset of new "features" is at z 2 = U^, which we require to be 
a large physical length: U^Roa/UkY^ 2 3> £ 2 S , which is equivalent to the original 
condition A5 3> R x given for the supergravity limit. Hence, we trust the supergravity 
approximation to describe the features near u — > Uk- 
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5 Analysis of equations: Stability 



A few words are in order to explain how we will address the issue of stability. We 
will not solve the eigenvalue problem, numerically or otherwise, to establish the 
four dimensional masses for fluctuations as being positive definite with a mass gap. 
Instead, we will simply show that to quadratic order in fluctuations, all actions are 
of the form (up to gauge) 



d 5 Z 



Fo&Qf) - ^1(6; Qf) ~(d t M) 2 + J2(d Xi Mf 



1=1 



-•F 2 (6 Q/)(%M) 2 - :F 3 (£; QfW 



(5.1) 



where M is a field describing the fluctuation, and 

Q/ = 0) >0. (5.2) 

This last statement is important, because it implies that anywhere the perturbative 
analysis is valid that (to linear order in Qf) 



Hi\ Qf) = Hi\ Qf = o) U + fr (d Qf r t ) \ Qf=0 j > o. 



(5.3) 



The inequality holds because if the term added to the 1 must be small. If it were not 
small, we would be forced to go beyond linear order in Qf, and hence the perturbative 
approach would no longer be trusted. Therefore it does not change the sign of any 
of the functions where the perturbative analysis is valid. Hence, the action that we 
have written had positive definite hamiltonian 



i=l 



+F 2 (Z;Q f )(d e M) 2 + Js(£; Q f )M 2 



(5.4) 



simply because it is a sum of squares. 

For this reason, what seems to be important is the presence of a cancelation 
between the DBI and CS action at the order \M\ (we expect such terms from the 
solutions above). 

Of course another interesting question is what happens outside of the regime 
of validity of the perturbative approach. This, however, is out of the scope of our 
present investigation, although we hope to address this in some future work. 
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5.1 DBI and CS equations of motion: x\ = (ttR x ) solution 

Of course to expand an action, we must expand about some solution to the equations 
of motion. For this reason, we briefly outline (and give more detail in appendix C) 
why £4 = (24 = irR x ) is still a solution to the equations of motion. One may 
address this simply by looking at the gauge invariant information in the equations 
of motion: namely the cusps. The cusps in (j)y,Ay,Gy and Cy cannot be removed 
with a coordinate transformations, as this would introduce delta functions into By. 
In appendix C, we show that the cusp in By, which is pure gauge, does not enter. 
We can read off the behavior around the cusps to be 



k = UM + n + 0(xl) 

1 {Qc9s) 2 

! V*ui(l-%)*\x A \ 
Fx = fi{u)--2 \- i + 0{x\) (5.5) 



2 (Qc9s 



2 



1 VSul (l - ^-V \x,\ 
F2 = h(u) + -4 \- -V + 0(xl) 

simply by comparing the delta function source terms and the coefficient (a function 
of u) of the d^ 4 term. The next order contributions are of order x\ because we require 
the functions be even about 24 = to obey the Z2 symmetry of the problem. 

Now we need to find the equations of motion for the embedding functions X M 
resulting from the action 



9sS_ 



B 



Qf 

In the appendix, we show that the only interesting equation of motion is the one for 



L at X4 = 0, similar conditions apply at X4 = ttR. x 
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X 4 , and we find that this becomes 



3 

U 4 



gsR D l 



«2 ( 1 



F4 



/ 



-±d. 



•1'4 



(#D4^ 

1 M4 - if 



V 



3 

# 3 



D4 



v- 



d X4 (|x 4 | ± \x A \) 



(5.7) 



where in the third line, we ignore the higher order in x^ corrections, as we will 
evaluate the derivative at x± = 0. Above we have also switched back to the more 
familiar g s Q c = 3i?|, 4 notation. 

The term with the ± comes from the CS term, and may vanish for the — sign 
choice above. We interpret this as putting a brane next to the backreacted branes. 
This tells that the equations of motion are satisfied for a brane placed directly on 
top of the other branes. If instead we had put an anti-brane, we would have found 
a constant force type potential, and we take that this is a solution too (although we 
expect an open string tachyon for small enough distances: our actions do not contain 
terms for strings ending on different branes). 

Further, we should note that the above is a gauge independent statement. The 
cusps in the functions Ai, C\, G\, are independent of the gauge choice 0, and 
these were the only functions that contribute above (see appendix C). Therefore, 
the leading \x±\ dependence is unaffected by small coordinate transformations. For 
this statement, it is important that B does not appear: it's values (cusps and all) 
are gauge dependent, while the other functions are determined in a gauge covariant 
way. Recall that while <f>i is gauge dependent, it's cusp behavior is not: one may not 
remove any part of the cusp without introducing unwanted delta functions into B\. 

This gives that to lowest order, original embedding solution is still a solution to 
the equations of motion. To truly consider the stability, however, we would like to 
know whether this extremum of the action is a maximum or a minimum. For this 
we will need to investigate the second order action about this point, and this will 



12 in any sense: small coordinate transformation or shifts of by a infinitely differentiablc 
globally exact form 
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involve second derivatives in £4, rather than just first derivatives. Hence the even 
functions that we were able to ignore in the above discussion will enter. 



5.2 Uk = decompactification limit 

We begin with the solutions for the decompactified case in the last section 

Fi(u,x 4 ) = uKi{q) q=— — -j 

{Qc9s) 2 

Mu,xA) = uK 3 (q), (5.8) 

where Ki are given in equation (14.51) . We will want to construct the second order 
action in x 4 and so we will need Ki to second order in x± oc q. We expand the Ki to 
obtain Fi and 0i and find 

, ^ 256 VSuhxA 768u 2 x 2 
F 1 [u, x 4 = uK x (q) = ——u - — — -j- + — — — 

1001 (Qc9s)^ l ^Qc9s 

F 2 m,x 4 = uK 2 {q) = -N 2 —u+ ' 4| +iV 2 — i 5.9 

9 (Qcgs)* 3 Qc9s 



' 4 2 (g^)! neg^ 

where we have dropped order 0(xl) and higher terms. 

We are now able examine stability of the decompactified limit by examining 
the second order action. We start by writing the pullback metric as a function of 

X A (x»,u) 

ds 2 p = e 2A (ri^dx^dx") + e 2G du 2 + e 2C d£l\ + e 2B (<^X 4 + d u X 4 du) 2 . (5. 10) 

We need the determinant of this metric to second order in X 4 , however it is easier 
not to expand the above metric completely, and simply realize that 

= (l + \<h ab ) + 0((X* f) (5.H) 

where g P Q is constructed by dropping the last term in ds 2 and h a b is the symmetric 
tensor defined by the last term in ds 2 . Evaluating this, we find 

- J d 9 ^^f p = (5.12) 
- J d ^e-^ A+4C+G (l + e^-^rT^XUX 4 + e 2 ( B ~ G ^d u X 4 d u X 4 
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again, dropping order (X 4 ) 3 and higher. We now take the expansion of the functions 
A = Aq + QfA\, ■ ■ ■ and evaluate to zeroth and first order in Qf, using the above 
expansions about X 4 = 0, and keeping only those terms second order in X 4 or lower. 



N 2 ^ + N 3 




1 






u 


20 J 




33^ 


i u 2 Xl 


80 , 


' Qc9s 



(5.13) 



45 2oy y 2 

To this we must add the term coming from the RR coupling. To do so, we recognize 
that 

A„ i]x4 = T^Fp—e 13 ^^ (5.14) 
9s 

satisfies the equations of motion for the nine form potential. We use the notation 
| £4 to mean that the index for X4 has been omitted. The factor of — g p in (15.141) is 
constructed using only the zeroth order in Qf metric pulled back. Further, this — g p 
only has corrections of order 0((X 4 ) 2 ), and so we may ignore them because of the 
|a?4| already multiplying yj—g p . Hence, to the order that we are working, 

- j 'A 9 = - J d^(T)u^\X 4 \. (5.15) 

Thus, for the correct orientation of the probe brane, the \X 4 \ term in the action 
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completely cancels 0, and the total action becomes 




9*u* 




1 5005 

Ar 8^ Ar 33 \ u 2 X\ 
^-jV + iVs— U-^ 5.16 
15 80 / Q c # s 



512 
5005 


+ iV 2 


9 


-N 3 ^ 
20 


Q/( 


512 
5005 


+ N 2 


4^ 
45 + 





l -d u x A d u x A 



Although we can at this point find equations of motion for the above action, and 
proceed with the analysis directly, we find it convenient to manipulate the above 
equation a bit more. For this, we note that we can redefine the u coordinate as well 
as the field X 4 . We find it convenient to do the following transformation 

u = u-\-Qf\\u 2 
du = (l + 2Qf\iu)du 

X\x^u) = (^+ l -Q f \ 2 ^jX\x^u) (5.17) 

(we do the u coordinate change first, and then the X transformation) and then for 
ease of notation we simply drop the A from the above. Again, we may only keep order 
QfOi lower in the above expansion. After doing this, we will introduce terms of the 
form Qf\ 2 f(u)X 4: d u X 4: = f X 2 f ' (u)d u ((X A ) 2 ) which we integrate by parts. This 



13 We interpret this as a brane next to the backreacted brane(s), rather than an anti brane next 
to the backreacted brane(s). 
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affects the coefficient of (X 4 ) 2 . After doing so, we find that the new action is 



/ 



„A ( Qc9s \ : 
- I d 9 £- [ 3 > 



g s u* 



512 AT 4^3 9 9, 

1 - Qf\ h N 2 - L - + N 3 Ai \u 

^ i 1 5005 45 3 20 2 1 ' 



1536 Ar 8^ AT 33 ll x \w 2 X 2 

+ - * + w »t + ^ - i Ai - A2 ) ») y^ xt 

The interpretation of Ai and A2 is that the correspond to coordinate transformations, 
and so are actually arbitrary and one may choose these. The numbers N\ and N 2 
are numbers that determine part of the profile of the backreaction of the branes, 
and so may be constrained by some physical boundary conditions. Here, however, 
we simply note that A2 may be chosen to eliminate the (X 4 ) 2 term completely. 
This then leaves terms of the form (1 — QfuCi) times terms present when Qf = 0. 
Therefore, we conclude that when the perturbative analysis is valid, all coefficients 
remain the same as the Qf = case. Because of this, one may simply argue that 
the hamiltonian of the above action is positive definite (it is a sum of squares times 
positive coefficients) for the range of validity of the perturbative analysis. Hence, we 
conclude that in the perturbative regime, the configuration is stable. 

This depended on the leading order cancelation between the DBI and coupling to 
the RR field. Other than this, the remaining terms were all present in the Qf — > 
limit (up to gauge). In such a case, all corrections that are order Qf cannot change 
the signs of coefficients, and so stability (in the range of validity for the perturbative 
approach) is preserved. We will see this again in the next section. 

As a curious note, with an appropriate choice of iVj and Aj, one can completely 
cancel the leading order in Qf contribution to the above action. 



5.3 Stability of the Sakai Sugimoto model 

Above, and in appendix C, we show that X 4 = 0, X 4 = ttR x is still a solution 
to the equations of motion resulting from the DBI+CS action. To evaluate the 
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second order action for fluctuations, we change to the radial coordinate z defined by 
v? = U K + Ukz 2 , and then to the "Cartesian" coordinates 

yi = z sin (jf^j V2 = z cos ( 5 - 19 ) 

which also allows for comparison with the analysis performed in [10]. Here the 
important point is that we chose the gauge B\ = G\, and so the only change in 
the z, £4 plane is by a conformal factor. Hence, much of the analysis of [10] follows 
through. We find that the metric in these coordinates is written 



ds 1 = e ZAl —— (rj^dx^dx") + e l ^R\^d£l\ (5.20) 

\-OD4 / 

3 

+e 2 Gl 4 SRdA 2 ^ _ h{z)y ^ dy 2 h{z)y 2^ dy 2 _ 2 h (z)y 1 y 2 dy 1 dy 2 ). 



9 V u 

where now all metric functions are written as functions of y± and y 2 , and we have 
defined the following functions 

h ( z ) = — 1 

z z \ u 

u = u{z) = (U K + U K z 2 )^ (5.21) 



z = z(yi,y 2 ) = \Jyl + yl- 



Here we have suppressed the factor of Qf for ease of notation, and will only rein- 
troduce it at the end. Taking the embedding yi(x /1 , it) one may compute the second 
order action the same way as the Uk = decompactification piece. One writes the 
line element as ds 2 = dsf + ds 2 where dsf is diagonal and ds 2 , is already order yf, 
and so again one finds that 

V^Fp = v^9p~o (l + + 0{{ yi f). (5.22) 

One may compute the DBI action easily now, 

-K 8 J d 9 ^^g~ p = (5.23) 

3 1 

3 g s J V 9 u 

lv? e , AM+G ^ {h{z){yl _ 2ymiJi) + y?) 
1 Uk 

14 



where we have defined d m y\ = y\. Further, the above function u(yi,y 2 ) still must 
be expanded in y\. The term V4 is the volume of the unit four sphere, and d 5 ^ = 
dtdx\dx2dx^dy2- 

We will integrate the term linear in y\ by parts, but first we find it convenient to 
introduce the following notation 



From the arguments in the last subsection, we expect the above combinations of 
fields to have the following behavior about yi — 



where in B we ignore higher corrections in y\ because its coefficient is already 0(yf). 
In the above, we have determined the expansion in A of order \y±\ by considering the 
argument in the u, £4 coordinates used to give equations (15.51) . and then changing 
coordinates to the 2/1,1/2 variables. 

In the following, we will have to evaluate u(yi,y 2 ) & t Vi = 0, and henceforth, we 
will call this function u y . Similarly we define h y = h(z)\ yi=0 . Plugging in the above 
to the second order action, and reintroducing Qf we find 




(5.24) 



(5.25) 



A = A (y 2 ) + -\ yi \ + A 2 (y 2 )y 2 1 + ■ ■ ■ 
B = 80(2/2) + ■•• 



(5.26) 




(5.27) 




9 u. 



2R 3 r 



(l + QfBoWd^y! 



1 u 

L UK 




where we define 



3 1 



2 K % RIJJIV A _ f 



(5.28) 



3 g s 
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as in the work of [10]. It is easy to read off the result in [10] in the Qf = limit; it 
is the top two lines of the right hand side. As in the last sections, we now add to 
this the contribution from J Ag. This is relatively easy to do, as we find 



-9 = ~RlMul ( 1 + + 0{y{ 

U K 



(5.29) 



and so to second order in y\ we can simply take 
A , -t^ 4 



T^Rhufrilvil (i + o{yl)) 

9s J 



(5.30) 



This exactly cancels the \y\\ term (for the — choice), as we have seen several times 
now (see appendix C for this occurring at the level of the equations of motion). 
Therefore, the full action reads 



K 8 j d 9 te-*^g~ p - K 8 Ja 9 



-T / d 5 £ 



{u 2 y + yl){l + Q f A ) 



2 R m { 

9 n„ 



l + QfBo)rrd^yid vyi (5.31) 



1 u 3 f 



1 u 

^j^h y y 2 d y2 A ) yl 



At this point it is sufficient to note that because all terms in the action were present 
before the perturbation, we expect that wherever the perturbative analysis is valid, 
the stability of the Sakai Sugimoto model is maintained. This is because whenever 
the perturbation is small, the above action yields a positive definite hamiltonian, 
and so all fluctuations will have positive energy. Further, as we have seen in the 
previous section, the equations admit a perturbative solution about u = Uk and so 
the perturbative analysis is valid from u = Uk up to u g-, where for the previous 
sections analysis gives a good approximation to the solutions. 



6 Discussion and outlook 

Here we will summarize our results. From the above calculations, we can see that 
when u is large enough, the solutions tend to that of the decompactified case. The 
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height of these functions all grow as u and so to stay in the perturbative regime, we 
require that 

1 4?r4 3 6Att 3 £ 3 s 

MQy<<1 ^„ << _ = _^ M 3 << _ 1 m 

There is a further requirement, that the supergravity approximation is valid. This 
gives a restriction 

u 3 V . 3 nNJ 3 



9 *Ud « 1 - i «^r- (6 - 2 » 

One may easily compare now and see which condition is more stringent, as all coef- 
ficients of g s and i s are the same. We find that generically 

1 R 3 n N C N 3 

< -> K -7^-4 (6-3) 
Q) g A s 64tt 2 v ; 

as we assume that N c is large and Nf ^ (see figure [7]). However, we note that for 
small g s the regime of validity of the perturbative backreaction (and so the validity 
of the probe approximation) becomes arbitrarily large. 




non perturbative Q f regime 



Figure 7: Range of validity for calculations 

In figure [7J we have indicated a possible M-theory (11D SUGRA) lift, although 
some caution is necessary, as no know lift D8 branes is understood in the context of 
11D SUGRA, at least for those described by the Romans type IIA. 

We are also left with some obvious open questions: 
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1. The topic of this paper has been the low temperature limit of the Sakai Sugi- 
moto model, and one may wish to know the qualitative differences between the 
low and high temperature limits. Further, one may hope that the analysis of 
the high temperature limit may be easier, as the D8 branes are transverse to a 
cylinder, rather than a cigar. 

2. It would be interesting to address the backreaction of flavor branes in other 
brane systems using the above techniques. While one may worry about the 
perturbation breaking down near the brane for codimension other than 1, one 
may trust the cancelation between the DBI and CS terms in the quadratic 
action for fluctuations. We believe this to be true for the following reason: 
for a section of brane near a smooth point in a manifold, it's backreaction 
(non perturbative contributions included) should behave just as the flat space 
case. In such a situation, a section of parallel probe brane near by feels no 
force on it because the charge and mass (per unit p volume) are the same. 
We may expect this to always be true. Further, in supersymmetric situations, 
the supersymmetry of the backreaction may be of some assistance in fixing all 
coefficients. We look forward to addressing these issues in some future work. 

Acknowledgments 

We are grateful for discussions with Oren Bergman, who was involved in the initial 
stages of this project. We would like to thank Ofer Aharony for discussions and 
comments on the draft of this paper. We also wish to thank Martha Merzig for help 
editing graphs, and Josh Davis for directing us to a useful reference. The work of 
B.B and J. S has been supported in part by the Israel Science Foundation, by a 
grant ( DIP H.52) of German Israel Project Cooperation D.I.P and by the European 
Network MRTN-CT-2004-512194. 

This material is based upon work supported by the National Science Foundation 
under Grant No. PHY-0455649. 



18 



A Massive type IIA equations of Motion 

We recall the following definitions 

F 2 = F 2 + MB 2 

F A =F A + -MB 2 AB 2 (A.l) 
F A =F 4 -A 1 AH 3 + -MB 2 A B 2 . 
The equations of motion for the action Sua m (12. 9ft are 

R, v + 2V At V,0 - ^H^ P°H vpa - e 2 ^^ f 2F„ <>F vp - l -g pv F 2 ■ PA (A.2) 
1 / ~ ~ 1 ~ ~ \ e 2 * 

R - 4< 90 . d(j) _ JL# 3 . H 3 + 4<T V^V^ = (A.3) 
V ai ( e - 2 *H aM + A aa F a ^^) - 2! 4! _ 1 4! e^fe-F (4) ...F (4) ... 



.MF^ - -MB aia i aia2ftft , - ■ r' Ii...li...F. u ... = (A.4) 

2 12 2-2!-4! v /= <? 



V Ql F a ^ - lF^ a2a ^4^ Q2Q3Q4 = o (A 5) 

Ql 31-4!^/=^ ^...^(4)... 

6 AAA '"BJ..I.. = (A.6) 



—B-F 2 - -B...B...F7 
2 4! 4 



--e~B...B...B...F {A) .., -M + *F 10 



2! -4! -4!v^ 

= (A.7) 

where again e takes values ±1. In the above, where we have written • • • there are 
indices contracted. To reintroduce the (sub) superscripts, one puts in a set of indices 
in the superscripts, and then puts the same indices in the same order in the subscripts. 
In the above, we note that M is constant, and may be considered piecewise constant 
in the presence of sources. 
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B Separating the equations 



Here we deal with the Einstein equations and the dilaton equation, and explain how 
to separate them. There are 5 Einstein equations, and one for the dilaton, and we 
name them 

EOM^ = 2d 2 J - (d u $f - 2d 2 u B - A (d u A) 2 - 2 (d u B) 2 - A (d u C) 2 
+2dJd u B - 2dJd u G + 2d u Bd u G 

+e (2G ~ 2B) ^2d 2 x J(u, xA) - (d Xi $) 2 - 4 (d X4 A) 2 - 2 (d X4 G) 2 - 4 {d XA Cf 

+2d x Jd X4 G - 2d x Jd Xi G + 2d X4 Bd X4 G - 2d 2 X4 Gj (B.l) 

2G _ 2C e 2A-B + 2G + 2C + ^ QfA 



-12e 



9s 



EOM tt = Ae^-^du (e^ G+B '^d u A^ 

+4e (2G-2B) eiB-G^Q^ ^-B+G-$)qA ( B .2) 



2A-B+2G+2C+U n A 
)2 4A+2G-4C+cj> o C V/^ 



— <5 e * — 2- 



9s 



+4e (^-2 B ) ( 4 (^ 4 A) 2 + 4(a c4 C) 2 + (d, 4 G) 2 + e^ 4 (e^fi^G - 0))) 

_ Q 2 e 4A + 2G-4C+0 + 2 £ (B3) 

EOM uu = A(A{d u A) 2 + A{d u C) 2 + {d u B) 2 + e^d u (e^d u {B 



p 2A-B+2G+2C+U n a 
. g 2 e 4A + 2G-4C+0 _ 2 £ ^//^ (B 4) 



<7s 
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EOM ss =-Ae( G - B+ ^d u (e^ G+B -^d u c) - Ae^ 2 ^ ^ B ~ G ^d X4 (e^ B+G -^d X4 C) 

+ 12e (2G-2C) _ g2 e 4A + 2G-4C+0 + ^_ V/ft ( B . 5 ) 

9 s 

EOM m = -Ad X4 Ad u A - Ad X4 Cd u C + d X4 dJ - d u Bd x J - d X4 GdJ (B.6) 

Each of the above equations is to be set to zero. The labeling we have used is that 
the subscript <p denotes to the equation of motion, tt denotes the time-time and 
XiXi Einstein equations (these are just one equation), xx denotes the x±x± equation of 
motion, uu the uu, ss the directions along the sphere, and m the mixed 11x4 equation. 
Further, we have taken A = A(x4) to be a function only of 24. 

We now wish to perturb the following equations about the background solution. 
For this purpose, we take the following expansion 



A(u,X4) = -ln\ —j— y- ) +Q f A 1 (u,x 4 ) 



B(u,x 4 ) = 7 ln I 7iwi J + ( 1 ~ ^f- ) + Qf B M O 




G(u,x 4 ) = -ri n (___J_)_li n (i_E|) + g /Gl(M7X4) 



C(u,x 4 ) = -7 In a a +ln(tt) +Q f C 1 (u,x 4 ) 



0(m,x 4 ) = - In ( I - 4 ln(u) + 2 \n(g s ) + Qffa (u, x 4 ) 

.Qc 3 9t 3 

It is now straightforward (and rather unilluminating) to expand the equations of 
motion and keep only the linear term in Qf. The only key point is that the source is 
already linear in Q f and so one plugs in the background fields only to the exponential 
e 2A-B+2G+2C+\<j> a pp ear i n g with A. Rather than writing this out explicitly, we will 
simply explain the steps involved needed to separate the equations. Henceforth when 
we write EOMi we mean the above equation of motion expanded to linear order in 
Qf. First, the most useful equation when expanded is equation EOM m flB.6p . and 
this becomes 

- -d X4 C x + l-d^G, - -d^A, + d X4 dJ L - jkii ( f 3 + ^|\ = 0. (B.8) 
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This may be integrated to give 



-Ci + ~G X - -A, + dj x - % ( " 3 + ^ + F(u) = 0. (B.9) 
u 2u u 4 u(u 3 - U K ) 



One may solve this for G\ and plug into the other equations. We will denote doing 
so as EC-Mild- One may easily solve for F(u) now, 

2EOM i \ Gl + £0M xa; | Gl = -4 V v ; v ; 



u(u 3 — U K) 



(B.10) 



and so 



F( "' = SR3T^)' (B11) 

However one can easily see that this perturbation is simply taking Uk — > Uk + &Uk 
and linearizing on 5Uk- This is because under this shift, neither 0i nor A\ nor 
C\ is changed, and so only G\ shifts in EOM m . The linear shift of G is given by 
3 u3-iri k ■ Thus, we may safely absorb F{u) into a shift into the definition of Uk- If 
need be, we may always reintroduce it by shifting equations that depend on B\ or 
G\ appropriately. Further, this is only a zero mode contribution (in £4) and so will 
leave unaffected much of our discussion. For these reasons, we take F[u) = for the 
time being, knowing how to reintroduce it later if need be. 

At this point we have eliminated 2 equations of motion at the cost of 1 function, 
which puts us on course to decouple the equations. 

Next, we make the simple observation that in all equations of motion EOM i \c 1 
only d u Bi and d\B\ appear. Thus, if we can solve for d u B\, we may eliminate Bi 
completely. We do so by taking EOM xx \ Gl — EOM uu \ Gl and solving this for d u B 1 . 
This combination still has a delta function, and so it is important at this step that 
A is a function only of £4 so that when the expression for d u B\ is substituted into 
d u d u B\ no derivatives of delta functions appear. 

We have now eliminated 3 of the 6 total equations, with the remaining combi- 
nations being EOM tt \ Glt9uBl , EOM xx \ Gl:9uBl = \EOM^\ G ^ 9uBl = EOM uu \ Gl:duBl , 
EOM ss \ Gi q uBi . However, here we find that 

-EOM tt \ Gl , 9uBl - lEOM xx \ GltduBl - ^-EOM ss \ Gu9uBl = (B.12) 
5 6 15 

Hence, we are left with only 2 independent equations for 3 unknown functions. This 
appears to be under constrained, however, these equations are actually equations 
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only of 2 linear combinations of the 3 functions. The decoupled combinations may 
be written 

-3EOM tt \ GuduBl - -EOM ss \ Gu8uBl = 



, 3( 4n 3 - f^)^ gsQcU^F\ 



5AuF 1 , 2u(Q c9s )IV3 



w 



Ki 

A (B.13) 



-£OM tt | GlABl - 3£OM ss | GlABl 
3^F 2 + 



, 3(Au 3 - U 3 K )d u F 2 , g s Q c u 3 d 2 X4 F 2 



u(u 



3 



4«(g c ^ s )3 V / 3 



2 



it 



(B.14) 



where 



A x = —F 1 + — F 2 -— 0! (B.15) 
5 10 1(T V ; 

C x = — F1 + -F2-— 0i. (B.16) 
1 10 5 10 r V ; 

We now turn to the question of fixing 0. For this purpose, we remember that we 
used the combination 



EOM xx - EOM, 



4 u^Q.gMji 2 (4gA ~ 5gugi - 5d u Bi + 12flA) 
3 fM 3 -f/3) 2 u 



' K 

4 



«(Qc9.)"V3 A (B17) 



to solve for d u B\. It is clear for our setup that <p must have some "kink" part in its 
solution to account for the delta function, as the only £4 derivatives that appear act 
on (p. However, by adding a zero —S/ud u (u J EOM m dx/±) to the above expression, 
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we find 



EOM xx - EOM uu - -d u (u I EOM m , dx A ) 

u 



4 u 3 Q c g s d 2 x ^ x 2 (u 3 - 7U 3 K ) djx 36uU 3 K 4> l 



3(« 3 -f/|) 2 u{u 3 -U 3 K ) (u 3 -U 3 K ) 



u{Q c g s )*V3 I u . 10 



T7 -r^^ + -(d u B 1 -d u G 1 ) (B.18) 



[u 3 - U 3 K ) 



(Qc9s 



Now it becomes clear how one may maintain continuity of the functions and at the 
same time separate the equations. We take B\ = G\ and then solve the remaining 
equation above. One may have guessed this gauge, as one can bring any two dimen- 
sional metric to a conformafly flat one. We do not impose this on the full metric, 
however, as the polynomials in u are easier to work with. 

One more comment is in order. If one wishes, one may linearize on a small change 
Uk — > Uk + SUk- Under this, B\ and G\ transform differently. This can give a new 
source term to the equation for 0i. However, this change only affects the zero mode 
(in x 4 ) of 0i, and hence will not affect the shape of 0i in the £4 direction. 



C EOM for DBI+CS: details of x 4 = solution 



J d 9 ^^i P + J A 9 , (C.l) 



Here we find the equations of motion for the embedding functions X M resulting from 
the action 

_ g s S B 
" Qf 

and explicitly show that £4 = (nR x ) is still a solution. We will find the equations 
of motion for the first part, and then turn our attention to the second part of the 
above action. First, we change frame by scaling the metric g s trin g = exp (§0)G to 
write the first part of the action 

S D = J d 9 ^^G~ p (C.2) 
The equations of motion for the fields X^(^) in the above action are 



5X^ 



-yf-GpG^u (v 2 p X u + Gf d ^ a d ^ b T u a ^j (C.3) 
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where objects with a p subscript are constructed using the pullback metric, and T is 
the full spacetime Christoffel connection. We wish to ask whether X % = £ l for i^x^ 
and X A = £4 =constant is a solution to the equations of motion. Consider first the 
X k : 

j& = -V=G;G u ^-j^d a G?y/=GfrX* (C.4) 

+Gp b d a X a dbX f3 -G lp {Gpa y p + G p /3 j0l — G a p^ 

where we use d a as shorthand for a partial derivative in £ a , and we have used the 
fact that our metric is diagonal. Choosing the X 1 = removes the derivative from 
the first part of the equation. Also, the fact that G is diagonal, and identical to the 
pullback metric for indices allows us to simplify the above further 

= —\/—GpG(p)ki ( / n diGp\/—Gp (C.5) 
\ V ~ g p 

+Gp b d a X a di ) X f3 -G' p t (5 ai Gu t N + 8/3iGu,M ~ G a p^) 

where % inside the parentheses are not summed. Note that the 5xi projects the last 
remaining metric down to the pullback metric, and we are left with 



which simplifies further to 



, / 1 .. .. 1 



(C.7) 



This is obviously zero: the metric is diagonal, so that 

- 1 v /"fti 1 riiiri _ n (r* c"\ 

p (j {p)ii — 1 ~ > Lr (p),i U (p)H < "p'Jfriiv — U \^-°) 

which then causes the first and third terms to cancel. 

Hence, we are left with evaluating the X A equation of motion: 

^G~ P G 44 (VjX 4 + G?d a X"d b XT%,) (C.9) 
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where we use the shorthand 4 to mean the X4 components. 

The first term vanishes as X4 =constant. The second term we evaluate similarly 
to the last discussion 

— a/ - G p Ga (VpX 4 + Gp b 'd a X L di,X M T A LM ^ 

— —yf—GpGu ^Gp^daX^dbX" -G 44 {S^G^ + 5 ui G^^ — G^ j4 )^ . (CIO) 

This time, however, the 5 ma and 5na give zero (as the X 4 =constant), and hence 
only the last term remains 

— V - 4 P (p) ab < 4 

^G P ^L=d Xi ^G~ p (C.ll) 

-Gr, 



= ^/^G p d X4 (AA G + G G + AC G ) 

where the subscripts are to denote the new frame that we switched to. Switching 
back to the string frame metric, we find 

yfT P e-+d Xi (~4> + G + 2A + 2Cj . (C.12) 

At this point we stop this analysis because we must also look at the equations of 
motion coming from the CS action, as this has a cusp as well. 

First, for the embeddings that we have chosen, the background form field is the 
following 



? 1Q = *M = *9l A 9 = n 9 ^e B \x 4 \ 



9s 9s 



-Qf 



A)12356789 — \J~9p e 1^1 (C.13) 

9s 

where £4 G {— irR x ■ ■ ■ ttR x }, producing both the positive and negative delta function. 
In the above, we must use only zeroth order metric functions, as the above statement 
is already linear in Qf. The function g p is the determinant of the metric setting 
X4 =constant, i.e. the pullback metric on the D8s. 
Next, we wish to consider the action 

S a = Ja 9 = ^J d^A Ml ... M9 d ai X kh ■ ■ ■ 0„„X " (C.14) 
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where the epsilon takes values ±1. This is very easy to vary w.r.t. the fields X : 



SS A 



+ ^dj(A Ml ... M9 )d ai X Ah ■ ■ ■ d ag X M *e a 



X 1 



(C.15) 



Taking / = i to be one of the directions along the world volume coordinates, we see 
that the first term and second term are identical. This is because in the first equation 
i and a\ must agree to give a non zero answer (when contracting the epsilon). Of 
course in the sum there are 8! occurences of this. Hence, the 9! cancels, and one 
simply gets diA ^ 9 \ X4 , where we use \x4 to indicate that the index for x 4 has been 
omitted. The second term is also equal to this, as the contraction yields a (9!). Hence 
the full equations of motion for the X 1 fields along the volume coordinates vanish. 

Next, taking the above equation of motion for the X4 direction, one finds only a 
contribution from the second part, i.e. 



Therefore, the full equation of motion for the field X 4 assuming that it is constant 
reads 



where we have restored a — sign earlier omitted in front of the DBI action. Also, the 
± is to be read whether we are putting D8 or D8 branes in the background. Both 
are linear order in Qf because to zeroth order none of the metric components depend 
on 24, and Ag is already linear in Qf. Hence, all other metric components are set 
to being their background values (except those with the £4 derivative). Considering 




(C.16) 



^g~ P e~*d Xi (--4> + G + 2A + 2C\ ± d^Ao-9^) 




(C.17) 
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the cusp solution near x 4 = one reads (factoring out Qf) 



-\f-9v 




d Xi (\x A \ ± \x 4 \) 



(C.18) 



where in the third line, we ignore the higher order in x 4 corrections, as we will 
evaluate the derivative at x 4 = 0. Above we have also switched back to the more 
familiar g s Q c = 3i?£) 4 notation. 
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